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Open Orbits in Satellite Dynamics 

F.P.J. Rimrott, F.A. Salustri 


Recent investigations into the effects of a constant atmospheric drag force on an essentially circular satellite 
orbit have led to an analytical expression for the resulting orbit, which is obviously not exactly circular, but so 
close to it that profitable simplifications can be made in its derivation. In the present paper several of the 
resulting type orbits are analyzed and compared. 


1 Closed and Open Satellite Orbits 

For a point satellite of constant mass m in the gravitational field p of a point master, the Kepler force (Chobotov, 
1996) is of magnitude 



0 ) 


If in additon to the Kepler force (1) there acts a negative transverse force of magnitude D, which we will assume 
to be constant, then Newton’s second law requires that 

-Pm _ m ^_ r 02j ^2a) 

r L 

-D = m(r'e + 2rQ) (2b) 

Depending on the magnitude of the term r in equation (2a) we may obtain for the orbit shape an ellipse, a circle, 
or a spiral. We shall refer to ellipses and circles as closed orbits, and to spirals as open orbits. 

The solution of equations (2) is exact if 

D = 0 and r & 0 (3) 

and results in an ellipse as orbit. Conditions (3) by the way, may also result in a parabola or a hyperbola; but 
strictly speaking, they are not orbits (Latin: orbis = ring). 

The solution is also exact if 


D = 0 and r = 0 and r = 0 (4) 

and results in a circle as orbit (Figure 1). 

The solution obtainable if 


D & 0 and 0 < mr « K (5) 

is approximate and results in a spiral orbit (Figure 2). 

The present paper is devoted to an investigation of the spiral orbits which result from conditions (5). 
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Figure I. A Circular Closed Orbit 



Figure 2. A Spiral Open Orbit 


2 The Ward Spiral 

By assuming that r can be neglected in equation (2a), equations (2a) and (2b) can be combined into 
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Equation (6) can be integrated in closed form. With the limits r=r and r = r t> for t = t and t = 0, the Ward spiral 
(Ward, 2000) results (line 1, Table 1) 
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This parameter equation can be changed into polar form (Figure 2) by using another approximation, viz. 


r o Q = 

as a consequence of which the Ward spiral takes on the form (line 1, Table 2) 
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A point satellite on a Ward spiral has an orbital energy, with eccentricity £ = 0 and semi-major axis a ~ r , of 
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It changes by approximately 
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per orbit, i.e., it decreases as the orbit spirals inward, with 
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It was Dr. C.A. Ward’s unique contribution to recognize the significance of the insignificance of the term r in 
equation (2a), which indeed is often so small that it can be successfully neglected. This situation is e.g. typical 
for satellites at high altitude on essentially circular orbits that are subjected to a constant atmospheric drag, 
which causes only minor orbital changes, i.e., such satellites remain on essentially circular orbits (Rimrott and 
Salustri, 2001). Ward’s approach leads to an integrablc expression, which eventually results in the Ward spiral 
(7). In the derivation process for this spiral there are numerous approximations, and the reader is well advised to 
keep this in mind. 

These approximations may be of a cumulative nature and lead eventually to a result that represents a poor 
approximation over all, or they may be of a corrective nature leading to a relatively good approximation 
eventually. 

We can change equation (6) into polar form if we approximate the angular speed either by 


0 = 


(13) 


or if we allow for the fact that the radius vector changes by only a small amount, by the somewhat simpler 
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(14) 


With equation (13) we obtain instead of equation (6) 
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Integrated between the limits r = r and 


r 


= r 0 for 0 = 0 and 0 = 0 we have 
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which appears in line 3 of Table 2. 

If on the other hand we employ equation (14) then equation (6) becomes 

dr_ = 2 Dr] n ~ ; 3/2 

dQ \i.m 

which results in 
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the same as equation (9). Equation (18) appears in line 1 of Table 2. 

Equation (6) can be differentiated with respect to time and then used to supply us with an idea as to the 
magnitude of mr . We have then 


m r = 6 —D 
K 


(19) 


an amount which fulfils condition (5) if the drag force D is small enough. 

Equations (16) and (18) are a first indication that the Ward spiral (7) may have competition from other spirals, as 
we will see in the subsequent section. 


3 Other Spirals 

The whole concept of open orbits is based on the assumption that r is so small that it can be neglected. This 
concept also obviously includes the case r = 0 . This in turn means that r = constant, which we can bring about 
by letting 
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instead of equation (6). Integration of equation (20) leads to 
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which is entered in line 2 of Table 1. Using equation (8) we obtain the polar form 
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which appears in line 2 of Table 2. 
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With the help of equation (13) we have 
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which when integrated supplies us with 
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which is entered in line 1 of Table 2. 

If on the other hand we choose to make use of equation (14) then 
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Equation (26) is entered in line 2 of Table 2. 

Of the many other possibilities to approximate the radius change rate let us begin with 
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Upon integration we obtain the spiral 
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and for the force 
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Equations (27), (28) and (29) have been entered in line 3 of Table 1. If we change the independent variable we 
can write 


dr _ 2D 3 
dQ p m 

and integrate to get 
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Equations (30) and (31) appear as line 3 in Table 2. 
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On the other hand if we begin with 
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we get 
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These three expressions are entered in line 4 in Table 
If we change the independent variable and write 
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we obtain 
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a hyperbolic spiral. Equations (35) and (36) have been entered in line 4 of Table 2. 
Lastly, let us use the approximation 
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which leads to 
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Equations (37), (38) and (39) appear in line 5 of Table 1. If we write, instead of equation (37), 
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and integrate we get 
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a logarithmic spiral. Equations (40) and (41) are entered in line 5 of Table 2. 
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Table 2. Polar Equations 
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4 Assessment 


Tables 1 and 2 have each a column 4, where the series expansions of the five different spirals are given. These 
series expansions turn out to be the same for all spirals, an indication that they are essentially equivalent for 
small values of the drag force D. Thus at first sight all five orbit spirals appear of equal suitability. Amongst 
them the Archimedean spiral appeals because of its simplicity. However, there is one property of the Ward spiral 
that makes it stand out from all the others, in that the Ward spiral predicts the applied torque correctly, as is 
shown in the following. The others do not, as the reader is invited to confirm! 

The satellite’s angular momentum, with £ = 0 and the Ward spiral’s equation (7), is 


H = m-J\Xr = m-J]Xr 0 - 
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Its derivative with respect to time is the torque 
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q.e.d. 

5 Conclusion 

The insignificance of the term mr in Newton’s second law in certain circumstances can be successfully exploited 
to obtain an analytic expression for open satellite orbits. In the present paper it has been shown that the spiral 
obtained by Ward shares characteristics with many other spirals, all of which lead to the same expression, when 
expanded into a series and broken off after the second term. As to the question which spiral represents an open 
satellite orbit the best, one might be inclined to give the nod to the Archimedean spiral because of its simplicity, 
until one realizes that the Ward spiral stands out from all others by a significant property. 
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